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We consider a chirally invariant lattice Higgs- Yukawa model based on the Neuberger overlap 
operator X^ "' . As a first step towards the eventual determination of Higgs mass bounds we study 
the phase diagram of the model analytically in the large iV/-limit. We present an expression for 
the effective potential at tree-level in the regime of small Yukawa and quartic coupling constants 
and determine the order of the phase transitions. In the case of strong Yukawa couplings the model 
effectively becomes an 0(4)-symmetric non-linear <r-model for all values of the quartic coupling 
constant. This leads to the existence of a symmetric phase also in the regime of large values of 
the Yukawa coupling constant. On finite and small lattices, however, strong finite volume effects 
prevent the expectation value of the Higgs field from vanishing thus obscuring the existence of the 
symmetric phase at strong Yukawa couplings. 
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I. INTRODUCTION 

Non-perturbative investigations of lattice regularized Higgs- Yukawa models as a limit of the elec- 
troweak sector of the Standard Model have been subject of many investigations in the early 1990's, 
see e.g. the review articles of Refs. [H, H, H, 0, H, S 13 • These lattice studies were motivated by the 
interest in a better understanding of the fermion mass generation via the Higgs mechanism on a non- 
perturbative level. In addition, the focus has been on the determination of bounds on the Higgs mass 
and the Yukawa couplings which translate directly into bounds on the - at that time not yet discovered 
- top quark mass. However, these investigations were blocked, since the influence of unwanted fermion 
doublers could not successfully be suppressed. Moreover, the lattice models of these studies suffered 
from the lack of chiral symmetry. The latter, however, would be indispensable for a consistent lattice 
regularization of chiral gauge theories such as the Standard Model of electroweak interactions. 

Here, we want to extend these earlier investigations in a new direction in order to overcome the 
previously encountered drawbacks by following the proposition of Luscher [8( for a chirally invariant 
lattice Higgs- Yukawa model based on the Neuberger overlap operator Q . Within this model an exact 
lattice chiral symmetry can be established while suppressin g th e fermion doublers at the same time. 
This is possible despite of the Nielsen-Ninomiya theorem [lOj, since the established lattice chiral 
symmetry is not the continuum chiral symmetry itself, but recovers the latter symmetry only in the 
continuum limit. We consider here a Higgs- Yukawa model including only the two heaviest fermions, 
i.e. the top-bottom doublet, and the Higgs fields. This simplification is reasonable, since the fermion- 
Higgs coupling is proportional to the fermion mass and hence small for the light doublets. We also 
neglect any gauge fields within this model, since they can be taken into account via perturbation 
theory. 

As a first step towards a numerical investigation of this Higgs- Yukawa model we begin by studying 
its phase structure. Here we present an analytical investigation of the phase diagram in the large 
iV/-limit following Refs. [ill. H^]. We refer the reader to these references for earlier works on lattice 
Higgs- Yukawa models. (See also Ref. [l3| for a first account of our work.) In the present paper we 
access the phase structure of the model at small and at large values of the Yukawa coupling constant, 
putting particular emphasis on the existence of a symmetric phase also in the strong Yukawa coupling 
regime. The latter strong coupling regime of a closely related, chirally invariant Higgs- Yukawa model 
in two dimensions was also studied in the recent work [l4j and corresponding Monte-Carlo simulations, 
performed in that model, support the existence of such phase [l5j . Extensions of our present paper, in 
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particular concerning the verification of the analytically obtained phase structure by explicit numerical 
simulations and addressing the question of lower and upper bounds on the Higgs boson mass, will be 
discussed in forthcoming publications. 

The outline of this paper is as follows: In Section [IT] we briefly describe the Higgs- Yukawa model 
considered here. In the following Section [Till we derive an expression for the effective potential in terms 
of the amplitudes of the constant and staggered modes of the Higgs field, which is a reasonable approx- 
imation at small values of the Yukawa and quartic coupling constants. We then present the resulting 
phase diagram in the large iVf-limit and determine the order of the occurring phase transitions. The 
phase structure in the regime of large values of the Yukawa coupling constant and arbitrary quartic 
coupling constants is then accessed by means of a different large iV/-limit presented in Section ITVl 
We show that a symmetric phase also exists at large Yukawa coupling constants. On small lattices, 
however, this symmetric phase is shadowed by finite volume effects preventing the expectation value 
of the Higgs field from vanishing. We then end with a short summary and outlook. 



II. THE MODEL 



Aspiring to investigate the Higgs Sector of the Standard Model of electroweak interactions, we 
consider here a four-dimensional, chirally invariant lattice Higgs- Yukawa model containing one four- 
component, real Higgs field $ and a number of Nf fermion doublets. The latter are represented by 
eight-component spinors tp^ , ijfo' with i = 1, ...,Nf. However, these Nf doublets are all degenerated 
within this model and correspond to the heaviest fermion doublet only, i.e. to the top-bottom doublet. 
This is a reasonable simplification due to the fermion-Higgs coupling being proportional to the fermion 
mass. We have introduced the fermion doublet number Nf nonetheless, because it will be possible to 
access the model analytically in the limit of large numbers of (degenerated) fermion doublets. 

Furthermore, there are also Nf auxiliary fermionic doublets X^ present in the model, which 
serve as a construction tool for the creation of a chirally invariant Yukawa interaction term. However, 
once the chiral invariance is established, these unphysical fields can be integrated out leading to a more 
complicated model depending then only on the Higgs field <E> and the Nf physical fermion doublets 
^W. The partition function of the given model can be written as 

r Nf r i 
Z= I D ^ {i) D ^ (i) D X (€) D X (€) exp (-£« - Sf n - S Y ) (1) 

•* i=i 

with the total action being decomposed into the Higgs action S$, the kinetic fermion action Sp, 
and the Yukawa coupling term Sy ■ It should be stressed once again that no gauge fields are included 
within this model. 

The four-dimensional space-time lattice, that the model is discretized upon, is assumed to have 
L lattice sites per dimension such that its total volume is V = L A . Here we allow for both, finite 
size lattices with even L £ N as well as lattices with infinite extension, i.e. L = oo, and we set the 
lattice spacing a to unity for convenience. The kinetic fermion action describing the propagation of 
the physical fermion fields ip^ is then given in the usual manner according to 

N t 

S kin =J2J2 #^V$ " 2px£> l n , m X% (2) 
i—1 n,m 

where the four-dimensional coordinates n, m as well as all field variables and coupling constants 
are given in dimensionless lattice units throughout this paper. Here, the (doublet) Dirac operator 
j)(ov) _ f)(ov) ^qjiov) j g gj ven by Neuberger overlap operator T>( ov \ which is related to the Wilson 
operator = i (v / + V 6) _ r V £v£ by 

V {o ^ = p 1 1 + -^L= I , A = - p, 1 < p < 2r (3) 
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with V^, Vj^ denoting the forward and backward difference quotients, respectively. In absence of gauge 
fields the eigenvectors and eigenvalues of the Neuberger operator are explicitly known. In momentum 
space with the allowed four-component momenta 



(-7T,7r]® 4 

{2-ku/L : n€N ,n< L}® 4 



p £ V = 

the eigenvectors of the doublet operator 2?(° 1 ') are given as 



u tk (p) 
Cu ek (p) 



for L = oo 
for L < oo 



C = ±1, e = ±l, k £ {1,2} 



(4) 



with u ek (p) denoting the usual four-component spinor structure 



« £fc (p) 



/f (ejla) fOT ^° and uek ^ = \H 



for p = 0, 



(5) 



(6) 



respectively. Here € (D 2 are two orthonormal vectors and the four component quaternionic vectors 
6, 8 are defined as 9 = (1, — if) and 8 = (1, +zf) = 8^ with f denoting the vector of Pauli matrices. 
It is well known that the eigenvalues v^{p) of T>^ ov ^ with Iin^^p)] ^ form a circle in the complex 
plane, the radius of which is given by the parameter p. These eigenvalues are explicitly given by 



v\p) 



P + P- 



ei^j p 2 + 2rp — p 
\/p 2 + (2rp 2 — p) 1 



pf, = siu(p M ), Pf, 



(7) 



The auxiliary fields on the other hand do not propagate at all and their contribution to Sp m is 
chosen such that the model will obey an exact lattice chiral symmetry. 

The Higgs field couples to the fermions according to the Yukawa coupling term 



N f 

n,m z—1 



(i - 75) 



(1 + 75) , t 



(8) 



where yjv denotes the Yukawa coupling constant and B n . m will be referred to as Yukawa coupling 
matrix. Here the Higgs field <E>„ is rewritten as a quaternionic, 2x2 matrix <j) n = $° 1 — i& n Tj acting 
on the flavor index of the fermionic doublets. Due to the chiral character of this model, left- and 
right-handed fermions couple differently to the Higgs field, as can be seen from the appearance of the 
projectors (1 ± 7s)/2 in the Yukawa term. Multiplying out the involved Gamma- and Pauli-matrices 
one can rewrite the Yukawa coupling matrix in the compactified form 



B m , n = 5 mtn ■ £($„), #($ n ) 



-$ 2 7 5 -M$i75 $° 1-^75 



(9) 



being block diagonal in position space. The model then obeys an exact, but lattice modified, chiral 
symmetry according to 



W _ 



Stp 



Sip® 



le 



75 ( 1 -j^' D(0V) ) V' W +75X W 



. S X® = ^75^ (<M, V (i) , 



2ie0 (10) 
= -2^ (11) 



with e denoting here the infinitesimal parameter of the chiral transformation. Since the (here omitted) 
lattice spacing a appears in front of the Dirac operators, this exact symmetry recovers the continuum 
chiral symmetry (after having integrated out the auxiliary fermion fields) in the continuum limit Q. 
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Finally, we use a slightly unusual notation for the Higgs action S$ given by 

S* = ~ KN ®l I^+A + *»-A] + E + A W ]T H^n ~ Nff (12) 



71. fl 



where kjv denotes the hopping parameter and Xn is the quartic coupling constant. This notation 
with Nf appearing in the quartic coupling term (which turns out to be more convenient for the later 
analytical considerations) can easily be shown to be equivalent to the more commonly used lattice 
version of the </? 4 -action 



n,fj> n n 

by rescaling the coupling constants A, «, y and the Higgs field ip according to 

,t. n \ A k y 

$ = C-(/3, Aat^^j, KA r = — , y N = —, 

where the constant C has to obey the condition 

C 2 - 2X N N f C 2 = 1 - 2A. 



(13) 



(14) 



(15) 



Furthermore, this latter action in Eq. Q13p can also easily be connected to the usual continuum notation 



f 1 t 1 



(16) 



which explicitly involves a bare mass mo and the forward difference quotient V/. This connection is 
established by scaling the field and coupling constants according to 



= 7-: A = C ■ A , «=— , y = y -C 



where the constant C has to obey the relation 



1 = C • 



2A C" 



(17) 



(18) 



For the further analytical treatment of this model we integrate out the fermionic degrees of freedom 
leading to an effective Higgs model given by 



Z = Jd$ exp(-S e// [$]) 



with the effective action S e f / [$] defined as 

Nf 



exp (-5* - - S> 



(19) 



(20) 



By applying some adequate substitutions the Grassmann integrations can be performed allowing to 
write the effective action S e f / [<&] in terms of fermionic determinants according to 



S ef f[&] = S*[$\ - N f • log [det (y N BV^ - 2pD^ - 2py N B 



(21) 
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III. LARGE JV/-LIMIT FOR SMALL YUKAWA COUPLING PARAMETERS 



In this section we will derive the phase structure of the introduced Higgs- Yukawa model in the large 
7Vf -limit for small values of the Yukawa and quartic coupling constants. The idea is to factorize the 
number of involved fermion doublets Nf out of the effective action S e f / [$] , since the integral over 
all Higgs field configurations in Eq. Q19p can then be reduced to the sum over all absolute minima 
of the effective action when sending Nf to infinity. This factorization can be achieved by scaling the 
coupling constants and the Higgs field itself according to 



VN 



!JN 



A/v — 

AT 



N 



k n = k N , $„ = 



(22) 



where the quantities yjv, \n, kn, and are kept constant in the limit Nf — * oo. 

One is thus left with the problem of finding the absolute minima of S e f / [<&] in terms of the latter 
quantities. In general the operators B and 2?( 011 ) do not share a common eigenvector basis making 
the analytical evaluation of the determinant in Eq. (f2Tj) impossible for general, space-time dependent 
Higgs fields. However, for sufficiently small values of the Yukawa and quartic coupling constants the 
kinetic term of the Higgs action becomes dominant allowing to restrict the search for the absolute 
minima of S e f / [<E>] to the ansatz 



$ n = $ • y/Nf ■ to + s ■ (-1) " 



(23) 



taking only a constant and a staggered mode of the Higgs field into account. Here $ € IR denotes 
a constant 4-dimensional unit vector (|$| = 1), and we will refer to m, s £ IR in the following as 
magnetization and staggered magnetization, respectively. 

For the actual evaluation of the effective action we use the fact that the matrix B now has a 
diagonal-plus-subdiagonal-block-structure in momentum space due to the chosen ansatz for the Higgs 
field according to 



v {ov) _ ^_ B ^ V (ov) _ 2p 



2p 



S(pi, P2 ) ■ S (P2) (<I) • (V^(p 2 ) - 2p) (24) 



+ a ■ 5( Pl , pa) • U( Pl , P2 ) ■ B<*) . (l)H( R ) _ 2 p 
+ 6{ Pl , P2 )-V^\ P2 ), 

where the diagonal part is caused by the constant mode of the Higgs field, while the sub-diagonal 
contribution is created by the staggered mode. In Eq. ([24]) this is expressed by p 2 denoting the shifted 
momenta p 2 = p 2 + (w, tt, tt, tt), where adequate modulo-operations are implicit to guarantee that 
p2 G T 5 - The matrices U(pi,p 2 ), T>( ov \p), and B^ are 8 x 8-matrices with the indices (,x€xk\, (, 2 e 2 k 2 
and denote the spinor basis transformation matrix 



U(pi, P 2) Cieiku(2e2k2 = [u Ciei * 1 (pi)] t U^ e! ' fca (p a ), 

the Dirac matrix 

and the Yukawa coupling matrix 



(25) 
(26) 



[u Cl£lfcl (p)] t B($)u C2e2fc2 (p) 



(27) 



= S, 



fcl ,/C2 



^!,e 2 %,c 2 -® + 5 eu - t2 {i( 2 5 (lX2 $ 1 +S Cu 



G 



respectively. Due to this diagonal-subdiagonal-block-structure the determinant in Eq. (j2"Tj) can thus 
be factorized by merging the four 8x8 blocks, which correspond to the momentum indices (p,p), 
(p,p), {p, p), and (p, p). Up to some constant terms, which are independent of <J>, we can thus rewrite 
the effective action as 



S eff [$] = S*[$]-JV>-log 

= &,>[$]- JV>- log 



det (p (ou) - |^ • B ■ (v {ov) - 2p 



0<p 3 < 



JJ det ( D<<"0(p) ® 0(°")(p) - ^M(p) 



(28) 
(29) 



where the restriction < p% < ir has just been introduced to prevent the double counting that would 
occur if one would have performed the product over all p G V after having merged the blocks. Here 
M.{p) denotes these merged, momentum dependent 16 x 16 matrices given by 



M{p) - ( M 2 > x (p) M 2 > 2 {p) 



with 



(30) 

(31) 
(32) 

(33) 
(34) 

The expression in Eq. (|29p can be written more compactly, taking the fact into account that the 
matrices involved in that expression are diagonal with respect to the index k due to Eq. (|26|) , Eq. ([27]) 
and 



M hl {p) = m-B { ri{$>)-(v {ov \p)-2p}, 

M 1 - 2 {p) = s-U(p,p)-B^(^)-(v {ov Hp)-2p 

M 2A {p) = s-U(p,p) (x>< OT )(p)-2p) 

M 2 - 2 {p) = m-B^(^)-{v {ov \p)-2 p y 



U(j>, fOcieiki.Caeaka = ^1X2 1 ^ 



(35) 

Since one easily finds that the determinant in Eq. (|29[) is invariant under the permutation p <-> p, 
one can extend the product in that equation, which is performed only over one half of the whole 
momentum space, again to the full momentum space V by factorizing out the identity &k\M- One 
then obtains for the effective action 



S e/ /[$] - S*[Q]-N f -log 



II det [D {ov \p)®V^\p) - V -fM{p) 



(36) 



with the definitions 

V^(p) = 8 klM -V^(p), M(p)=5 klM -M(p), and M a >\p) = 6 klM ■ M a > b ( P ), (37) 

where a, b E {1, 2}. Selecting a special order for the indices (e according to {++, H — , — h, } the 

latter four 4x4 matrices are explicitly given by 



M 1 ' 1 {p) = m 



M'' 2 {p) 



( $°w+(p) i&u]-(p) (z$ 3 - &)uj-(p) \ 

(il> 3 - $ 2 )cj+(p) 

(i$ 3 + $ 2 )w-(p) $°w+(p) -id 1 ^-^) 

+ &)uj + (p) -i$ 1 w+(p) 4°cj-(p) / 

/ i$ 1 w+(p) $°w~(p) (i$ 3 - $ 2 )u+(p) \ 

$°w+(p) ^^-(p) (i$ 3 - ¥)uj-{p) 

{i¥ + $ 2 )w+(p) -i$ 1 w+(p) $°^~(p) 

V (i$ 3 + <l 2 )cj-(p) $°w + (p) / 



(38) 



(39) 
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where the abbreviation w e (p) = v e {p) — 2p was used. The remaining matrices M. 2,2 (p) and A4 2,1 (p) 
are obtained from M 1 ' 1 {p) 1 M 1 ' 2 (p) by interchanging p and p. Using some algebraic manipulation 
package, the determinant of the 8x8 matrix in Eq. (|36p can be computed leading to the final expression 
for the effective action 



pEV 



v + (p)\-\v + (p)\ + P(m 



(to 2 - s 2 ) ■ \v + {p) - 2p\ ■ \v + (p) - 2p\^j 



+ m2 %{ - M ■ \ v+ b)\ K(^) - M ■ \» + (p)\j 



(40) 



With the ansatz in Eq. (f2"3")) the Higgs field action S$ can also be written in terms of the quantities 
to and s. One easily finds 



5$ = N f ■ L 



TO 



(TO 2 - 5 2 ) 



+ m + s + Aat(to + s + 6m s 



,2„2 




(41) 



Two remarks are in order here for the orientation of the reader. 

(I) The resulting phase structure in the large iV^-limit can now be obtained by minimizing the 
effective action with respect to m and s. In principle one could derive the corresponding phase 
diagrams for all values of the quartic coupling constant Xn > 0. However, as one can easily find from 
Eq. (|15p the case An > 0.5 corresponds to the strong self-coupling regime A ^> 1 of the physically 
underlying </? 4 -theory given in Eq. (|13[) for large values of Nf. In that regime it is no longer reasonable 
to evaluate the effective action due to the strong self-interaction of the Higgs-field in that case. We 
therefore restrict the allowed range for the quartic coupling to < An < 0.5, which corresponds to 
the weak self-coupling regime of the physical model in Eq. (|13p . 

(II) The sum over all allowed momenta V in Eq. (|40p becomes a four-dimensional momentum 
integral over V for L = oo according to 



-y 

pEV 



pEV 



d 4 p 

(2^ 



(42) 



which was actually used in the numerical evaluation of the effective action. 



Ajv = 0.1 



Ajv = 0.3 




KN 




AFM 



VN 



2 4 6 8 

VN 



FIG. 1: Phase diagrams with respect to the Yukawa coupling constant j/jv and the hopping parameter ftjv for 
the constant quartic couplings An ~ 0.1 (left) and An = 0.3 (right). The black line indicates the first order 
phase transitions. Both phase diagrams were determined for L = oo. An explanation of the occurring phases 
is given in the text. 
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We now present the phase diagrams for Ajv = 0.1 and A at = 0.3 in Fig. [T] These phase diagrams 
were calculated for an infinite lattice, i.e. for L = oo. Here we distinguish between the following four 
phases: 

(I) The symmetric phase (SYM): m = 0, s = 
(II) The ferromagnetic phase (FM): m ^ 0, s = 
(771) The anti-ferromagnetic phase (AFM): m = 0, ,s ^ 
(7 V,) The ferrimagnetic phase (FI): m^0, s =^ 

In both cases, i.e. Xn — 0.1 and Ajy = 0.3, one finds a symmetric phase approximately centered around 
kn = at sufficiently small values of the Yukawa coupling constant y^r, as one would have expected, 
since the model becomes the pure </> 4 -theory in the limit —> 0. From the same consideration one 
would also expect the accompanying phase transitions to be of second order. This is indeed the case 
as can clearly be seen in Fig. [3 showing the expectation values of the amplitudes m and s for different 
values of y n as obtained in the minimization process. With increasing y^ the symmetric phase bends 
downwards to negative values of the hopping parameter kn, unless it either hits a first order phase 
transition to an anti- ferromagnetic phase (black line in Fig. [I]), the order of which can be determined 
from Fig. [3] (this is the case for Ajv = 0.1), or it eventually goes over into two FM-FI and FI-AFM 
second order phase transitions, which is the case for A at = 0.3. 

Here we present only the expectation values of m and s for Ajy = 0.3 and not for An = 0.1, since 
the latter plots would not provide qualitatively new information to the reader. 

Interestingly, the ferrimagnetic phase (FI) exists in both presented scenarios, i.e. for Aat = 0.1 and 
\n = 0.3, even deeply inside the anti-ferromagnetic phase region in the neighbourhood of the first 
order phase transition boundary. However, due to the small expectation value of the amplitude m of 
the constant mode (see Fig. [5]) it is questionable whether this ferrimagnetic phase will be observable 
in corresponding numerical simulations. 



IV. LARGE iVf-LIMIT FOR LARGE YUKAWA COUPLING PARAMETERS 

In this section we will examine the phase diagram of the considered Higgs- Yukawa model in the 
regime of large values of the Yukawa coupling constant j/at and for arbitrary values of the quartic 
coupling constant Aat > 0. This will be done in three steps. Firstly, the effective action is expanded in 
powers of the inverse coupling constant 1/yjv- Taking only the first non- vanishing contribution of this 
power series into account and performing the large Yy-limit in such a way, that the amplitude of the 
Higgs field is fixed, the model then effectively becomes an 0(4)-symmetric, non-linear sigma-model 
up to some finite volume terms. Finally, the phase diagram of the latter sigma-model is determined 
by an additional large TV-limit, where N denotes here the number of Higgs field components. 

For an evaluation of the effective action it is crucial to pay special attention to the fermion doubler 
modes 



which we will refer to as 7r-modes in the following. Given these 120 modes one can define the corre- 
sponding projection operator 



projecting to the sub-space — span(Q,r) spanned by Q^. Using this notation one can easily 
establish the very helpful relation 



Qn = {* p Cek : P„ e {0, tt}, p 0, C, e = ±1, k e {1, 2}} 



(43) 




(44) 



det (E (1 - P^) + P*FP V ) = dct ((1 - P w ) E(l- P n ) + P^FP^) 

= det ((1 - P„) E + P-nFP-n) = det' (E) ■ det* (P) (45) 
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KN 

FIG. 2: Expectation values for the amplitudes of the constant (m: black curve) and staggered (s: gray 
curve) modes for several selected values of the Yukawa coupling constant yjv and a constant quartic coupling 
A at = 0.3. The results were obtained for L = oo. 



where E and F are arbitrary operators defined on the same space V as 2?( ou ' and B. Here the 
expression det* (F) denotes the determinant of F with respect to the sub-space V v and det' (E) is the 
determinant of E with respect to the complementary space V/V n = span(Q/Q 7r ), where Q denotes 
the full set of all modes. Using Eq. (145)) several times one can rewrite the effective action according to 



elf 

e 



det (y N B (V ou > - 2p) - 2pV ( - ov ^ (46) 
= (-4p 2 ) 120 ■ det' (y N B' (v l{ov) - 2pl') - 2pV' {ov) ) 

= (-4p 2 ) 12 ° • det' (y N ) ■ det' {v l{ov) - 2pl') • dot' (b' - ^V' (ov) (v' {ov) - 2pl'J 

= Const • det [ B - (B - 1) P n - —A] 
\ Vn J 

= Const • det (B) ■ det (l - (l - B^ 1 ) Pj ■ det ( 1 - — B^A [l - (l - P -1 ) PJ" 

\ UN 

where T>' l ' ov \ B' , and 1' denote the restriction of the operators V^ ov \ B, and 1 to the sub-space 
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V/Vtt- This restriction is introduced, since it guarantees V l<yCV ^ — 2 pi' to be invertible. The operator 
A is then defined by extending the domain of the inverse 

oiV l(°v)_ 2pl , 

again to the full space V by 

inserting the projector 1 — P n according to 



A = V 



(ov) 



V 



/(ov) 



2pl' -(1-Pk) 



(47) 



which is well-defined and finite over the whole space V. The last determinant in Eq. ([46]) can further 
be reduced by using the result 

[i - (i - b- 1 ) p v y l = i-p^+p^ii-p^+p^p^y 1 p„ 

- (l-P w )B- 1 P 7r (l-P w + P w B- 1 P w )~ 1 P n (48) 
and by applying again relation (|45[) leading to the compact notation for the effective action 



S e ff[$\ = S* - N f -logdet (B) -N f -logdet* (B' 1 ) - N f -logdet* ( 1 + ^-F[$] 

\ UN 

- Nf ■ logdet (1 - —A ■ B- 1 

V VN 



(49) 
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with the abbreviation F[$] defined as the somewhat lengthy expression 



m = 



VN 



B^AB^Pk [l-P w + P„B- l P„] 1 



(50) 



However, the latter determinants det* only give rise to some finite volume effects, since these determi- 
nants are only performed over the 120-dimensional sub-space V n . Their contributions to the effective 
action do therefore not scale proportional to L A as the lattice size increases in contrast to all other 
appearing terms. We will come back to discussing these finite volume effects later. Here, we will first 
continue with the evaluation of the last term in Eq. (|49|) by rewriting the corresponding trace as a 
power series in the inverse coupling constant 1 jy-^ according to 



Trlog ( 1 - —A-B- 1 

VN 



- Tr E 



2 r f p 



[AB-r 



(51) 



and by eventually cutting off this power series after the first non- vanishing term, which is well-justified 
for sufficiently large yjy. For our purpose of establishing the desired connection to a sigma-model it 
is most convenient to evaluate these expressions in position space. Here the matrix B~ x is block 
diagonal and explicitly given by 



B- 1 = Bt . (£?£}t) \ B-* n = 8 m , n 



, 2 ) 

n\ /) 



(52) 



where the notation ($*)° = (&n) 1 = — was used and B was defined in Eq. ©. In position 
space the matrix AB~ X can hence be written as 



[ab ]„ iri2 

p£V C,ek 



e ip ni u Qek ( p j a e Qj) e -ijma [ u Ctk ^ j t 



(53) 



71 E E « e b)e lp(ni -" 2) ^ fe (p) (B^(KJ\<S>n 2 \ 2 )) , , , \u^' k \p) 

Li — — \ / Cek.C'e'k' L 



pev f 

C'e'k 



with B^ as defined in Eq. (|2"T|) . The scalars a e (p) denote the eigenvalues of the anti-hermitian 
operator A corresponding to its eigenvectors ^p^^ and are explicitly given by 



iJRB a e (p) 



: pe V, is c (p) = 2p 



(54) 



The result for the trace of the operator AB 1 is then directly found to be 

Tr[AB-^ = ij^^Tr 8x8 [|$„|- 2 ^)sW($; 
n peV 

which can be generalized to the trace of the r-th power of AB^ 1 yielding 



(55) 



Tr [AB- 1 ^ = Tr ' 



8x8 



JL p ipi{ni-n i+1 ) 
II £ 4 \*r 



J- 2 A( Pt )(B^(K, +1 )) U{ Pu p i+1 ) 



(56) 



where p n +i is identified with pi, and x n +i with x\, and the expression A(p) stands for the diagonal 
matrix 



(57) 



At this point we refer the interested reader to Appendix [5] for the details of this calculation in order 
to sustain the readability of this text. 
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However, it turns out that the evaluation of Eq. (|56[) becomes much easier, if one inserts the identity 
U(pi,0)U(0,pi) at some proper places. The remaining 8x8 trace can then be simplified to 



Tr 



8x8 



r 

1[A( Pi )(b^(^ +1 )) Ufa 



Pi+lj 



Tr : 



8x8 



n^ (o) fe)(^(c, +1 ) 



(58) 



where the representation of the Yukawa coupling matrix can directly be taken from Eq. (127p 

and A 1 - ^ (p) is given by 



A^( P ) = U(0,p)A( P )U(p,0) 

a+ {p) _ ( -Je \ ^ f Po -p& 



Vie 



(59) 



where the relation a + (p) — —a~(p) has implicitly been used. Due to the insertion of the spinor basis 
transformation matrices U(pi,0) and U(0,Pi) the sums over the momenta in Eq. (|56[) factorize now 
according to 



Tr [AB- 1 ]' 



E 



8x8 




oipi(ni-m+i) 
L4 



(60) 



where each momentum sum is a four-dimensional Fourier transform of an anti-symmetric and purely 
imaginary summand, hence yielding real values. With the definition 



pev 



L 4 



(61) 



the hermitian matrix T n , m appearing in Eq. (|60j) can compactly be written as 



Tn 



$° n fe + i $i n r -i^fe-^r (**«,- $£,)r H^ + ^)fe 



$ 2 )r n 



(62) 



with the abbreviations T M = r^(An) and An = n — m. Therefore, the first order summand of the 
power series in Eq. ([5T|) reading 



Tr [AB- 1 ] = Tr sx8 [T n ,„] = 



(63) 



is identical to zero and the first non-vanishing contribution is the second order term, which can be 
evaluated by explicitly computing the 8x8 trace, yielding 



Tr [AB = E Tr 8x8 K.i,"2^2,ni] 



E 



ni n 2 



\T(An) 



(64) 



Cutting off the power series in Eq. (|5ip after this first non-vanishing term, which is well justified for 
sufficiently large values of the Yukawa coupling constant, the effective action can be written as 



s eff m = s* - N r (x;mi*«i 8 ) 



E i r ( An )i' 



(65) 



N f • logdet* (B^ 1 ) - N f • logdet* 1 + -t-F[§>] 

Vn 



2p 
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where the matrix has been defined in Eq. (J5T 

Some remarks concerning the remaining determinants in the latter result are in order here for the 
orientation of the reader. Here det* denotes the determinant over the sub-space V w , which has the 
dimension 120. In contrast to all other terms appearing in the effective action these determinants are 
not proportional to L 4 . They are therefore suppressed as the lattice size L goes to infinity. Moreover, 
the very last term in Eq. f|65[) even vanishes on finite lattices when the Yukawa coupling constant yjv 
becomes large. This is in contrast to the determinant det*(-B _1 ) being independent of ijn- However, 
it is nevertheless quite instructive to consider these finite volume effects in more detail. This can at 
least be done for the first determinant det*(B^ 1 ), which can be exactly evaluated for the ansatz given 
in Eq. (|23|) taking only a constant and a staggered mode of the Higgs field into account. In that case 
the inverse of B can also be described in terms of a constant and a staggered mode according to 

_ 1 / ^2 n M \ 777 S 

$„/|$„| 2 =$-^V • m + s-(-l)" , m=— -, ~s=-„ - (66) 

\ / 7T7 Z — Q z _Q — 777 z 



which allows to determine the desired determinant in a similar manner as described in Section IIIII 
yielding 

log det* (B" 1 ) = -60 log (Nf) + 8 log \fh\ + 56 log |m 2 - s 2 1 . (67) 

The obvious asymmetry in m and s is caused by the fact that the 8 zero modes vp°^ efc are not 
included in the sub-space V„ . The effect of the latter terms and especially the asymmetry in m and s 
is clearly observed in corresponding Monte-Carlo simulations [l(| on small lattices and large values of 
the Yukawa coupling constant yjv- Moreover, the result in Eq. (|67p would also hinder the expectation 
value of the Higgs field from vanishing, thus obscuring the potential existence of symmetric phases at 
large yn on small lattices. However, as the lattice size increases these finite volume effects eventually 
disappear. In the following we will therefore neglect the det* terms in the effective action (|55|) . which 
is well justified on sufficiently large lattices. 

To establish the announced connection to a sigma-model we now consider the large iV^-limit where 
the coupling constants scale according to 

Vn = Vn, Vn = const, \ N = \ N = const, k n = k N = const, (68) 

N } N f 

and for the Higgs field we consider an ansatz in which the amplitude of the local vectors $„ is fixed 
to (f S 1R according to 

$« = \fNf~ ■ <p ■ a n , \a n \ = l (69) 

where the four-component, space-time position dependent unit vectors a n are arbitrary. In this setting 
the contributions to the (reduced) effective action are either of order O(Nf) or 0(1). Considering 
only the leading order terms, for which the fermion doublet number Nf can be completely factorized 
out, then allows to fix the Higgs field amplitude if by the determination equation 

= -4 ■ -L + 1 + 2A W • (ip 2 - 1) . (70) 

With this fixation of the Higgs field amplitude the model in Eq. (f6"5)) becomes effectively a non-local, 
four-dimensional, non-linear sigma-model in the large TVf-limit given by 

Seff = - £ <{f„ 2 "<^ 2 (71) 

711,712 

with the effective, non-local coupling matrix 



±4 



tL = §^ |r(An)| 2 + K N ■ V 2 ■ J2 ^»,v (72) 



VnV 
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Here the notation "non-local" simply refers to the fact, that the field a n at any lattice site n couples 
itself to any other site of the lattice. This leaves nevertheless open the possibility that the interaction is 
local in a field theoretical sense with exponentially decaying coupling strength [13] • We did, however, 
not investigate the question in this paper, since eventually we are mostly interested in the small 
Yukawa coupling region. 

Basically, the outcome in Eq. ([72]) reproduces the result which was found for a Higgs- Yukawa model 
based on Wilson fermions [HJ with the only difference that the coupling matrix in that case consisted 
only of nearest-neighbour couplings. 

The phase diagram of the obtained sigma-model (|71[) can be determined analytically by an additional 
large iV-limit with TV denoting here the number of components of the vectors a n . The first step towards 
this evaluation is to remove the restriction \a n \ = 1 by introducing an auxiliary one-component, real 
field A„. This can be done at least in two ways. One can either encode the restriction |er n | = 1 as a 
(^-function [l8| written in terms of an integration of exp(iA n (|cr n | 2 — 1)) over A„, or alternatively, one 
can address this restriction by introducing the field variables X n as Lagrange- multipliers [l9| . Here 
we follow the latter approach which leads us to the extended action 

a i = ^{eE -<U ■ < ■ < + E A » • (E [<] 2 - 1) } (73) 

the minima of which can now be searched for without having to consider any restriction on the Higgs 
field amplitude. Here, an additional parameter was introduced. For = 1 the given action 
corresponds to the prior form of the action. This new parameter is inserted, since it will allow to 
factorize a factor N out of the action as required by the large N approach. This can be achieved by 
scaling according to 

tjsr = -tt, t N = const, (74) 

where we choose the setting tjv = 4, since this recovers our actual effective sigma-model at N = 4. 

The remaining problem to solve is to find the minimum of the action S[a, A]. However, it is well 
known from investigations of pure sigma-models that the phase transitions of such models cannot 
be correctly determined by evaluating the effective action S[a, A] in Eq. (|75|) directly by restricting 
the consideration to only some selected modes of the fields a and A. (Doing so would yield a first 
order phase transition at kn — 0.) This is in contrast to the situation we discussed in Section IlIII 
Instead, we first integrate out all modes of all N components of the field a except for the constant 
and staggered modes. This can be done by taking only the constant mode of the auxiliary field X n 
into account, i.e. A„ = A = const. Doing so reduces the action S[a, A] to 

1 f N 

S[m\s\\] = -ln[det'(-^ + A)] _Ar/2 + ^<j ^ [m*] 2 • (0 | - K e// + A| 0) 



i=l 



N } 

EM •< 7r l- Ke// + A l 7r >- i4A ' 
%=\ J 



(75) 



depending only on the real scalar A and the amplitudes m 4 , s l of the constant and staggered modes, 
respectively. Here the notations |0) and \tt) were used, denoting the constant and staggered modes 
(normalized by a factor 1 / y/lA) according to 

\k) ^je*" (76) 

being eigenvectors of K e ^ and det' is the determinant neglecting the two latter modes. For conve- 
nience, the introduced short-hand notation = (0,0,0,0) and n = (n, 7r, 7r, tt) will also be applied in 
the following where it is unambiguous. 

One remark is in order here for the orientation of the reader. The performed Gauss-integrations 
are only well-defined, if the involved eigenvalues of the operator —n e ^ + A are positive, which is not 
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guaranteed at this point. However, this step will be justified (and made more precise) a posteriori 
when a certain value for the scalar A will be assumed by solving the resulting gap equations. Here we 
will first continue with this formal expression and postpone its further discussion to the end of this 
section. 

To evaluate this latter determinant, the eigenvalues of the matrix need to be known. The 
eigenvectors are simply plane waves with wave vectors k S V and one easily finds the corresponding 
eigenvalues according to 



J.kn-2 



4 

2R,Nip 2 cos(fc M ) + 



16p 2 



q(k) 



where q(k) denotes the eigenvalues of the matrix |r(An)| 2 given by 



p = k-p. 



(77) 



(78) 



For the numerical evaluation of this quantity it is useful to use some symmetries of q{k). One has 
q(k) = q(k') at least, if k' is a permutation of the components of fc, or if k' = for all /i. 

Now we can search for the absolute minima of the effective action in Eq. (|75p . For this purpose we 
relate the amplitudes m 1 , s l to the values of the overall magnetization m and staggered magnetization 
s, respectively, according to 



■ rn 



and 



(79) 



With this notation one directly obtains from the effective action in Eq. ([75]) the following expression 
in terms of the quantities m, s and A 



Sim, s, A] = — Tr' In \-n eff + X] + ^- ■ m 2 ■ L 
2 L 1 t N 



=5^2 g(7r) + A ) - —L X 



fe^ ) + A 



tN 



(80) 



where the summation over the coupling matrix components has been performed by using Eq. {77J with 
the settings k = (0, 0, 0, 0) = and k = (tt, it, it, tt) = it, respectively Analogously to det', Tr' denotes 
the trace neglecting the modes |0) and \tt). We can now derive the corresponding gap equations by 
differentiating with respect to m, s, and A leading to 



m ■ 



16p 2 



9(0) 



(81) 



2 16p' , v 

y N v 



(82) 



2 i 2 

m + s 



t_N_J_ 

4 L 4 



E 



-k N Lp 2 ^2 cos(fc M ) 



(83) 



Equation (|81| implies that m or the given argument within the square brackets has to vanish. An 
analogous observation can be drawn from Eq. (|82|) . For the investigation of the phase structure we now 
consider two different scenarios for the amplitudes m and s, namely a ferromagnetic phase (m =/= 0, 
s = 0) and an anti- ferromagnetic phase (m = 0, s ^ 0). For each of these cases we can then derive a 
self-consistency relation: 
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FIG. 4: Phase diagrams for L = oo with respect to the Yukawa coupling constant y./v and the hopping 
parameter Rn for several selected values of the quartic coupling constant Ajv. The presented phase structure 
was determined for e = lO^ 1 , while the black lines show the phase transition lines obtained for e = 10 -3 . An 
explanation of the e-dependence of the presented results is given in the text. 



1. For a ferromagnetic phase (FM) (m 7^ 0, s = 0) one obtains from (|8"Tj) 

A = 8h N <p 2 + ' 9(0) 

and hence the following self-consistency relation 

%N<P 2 E ^ ~ cos(fc M )) 



< m =1 ——r 7 



4 L 4 



^4 ( ? (0) - 



Wm(*0 

2. For an anti- ferromagnetic phase (AFM) (m = 0, s ^ 0) one obtains from 

A = oK^ip + ^ — t ' 9 (tt) 
and hence the self-consistency relation 



(84) 



(85) 



(86) 



-kat^ 2 V (1 + cos(fc M )) + (g (tt) - g(fc)) 



W s (fc) 



(87) 
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FIG. 5: Expectation values for the amplitudes of the constant (m: black curve) and staggered (s: gray curve) 
modes for several selected values of the Yukawa coupling constant j/jv and the quartic coupling parameters 
Ajv = 0.0 and Ajv = 0.1. The results were obtained for L = co. 



Three further remarks shall be given here. 

(I) The equations and (|57|) are denoted as self-consistency relations because the assumption of 
a (anti-)ferromagnetic phase becomes inconsistent, if the resulting value for m 2 (or s 2 , respectively) 
becomes non-positive. If both assumptions become inconsistent simultaneously, this corresponds to 
a symmetric phase (SYM) with m = s = 0, while the case m 2 > and s 2 > is denoted as a 
ferrimagnetic phase (FI). 

(II) For the ferromagnetic phase the choice of A according to Eq. (|84p justifies the integration 
performed in Eq. (|75[) a posteriori, because it sufficiently shifts the eigenvalues 2W m (fc) of the matrix 
— n e ff + A to make all of them positive, except for the constant mode (fc = 0) which was excluded 
from the Gauss-integration. 

(III) For the anti-ferromagnetic phase, in contrast, choosing A according to Eq. (|5rj|) does not 
guarantee all eigenvalues 2W s (fc) of —n e ^ + A to be positive. The Gauss-integration in Eq. (|75p can 
therefore only be performed for all those modes ^ k ^ it which fulfill W s (k) > e with an arbitrary 
lower bound e > 0. The details of this statement are presented in Appendix [B] The results of this 
more careful consideration are already presented in Eq. (|55|) and Eq. ([57)) . The only difference to the 
naive result is that the set over which the sum has to be performed is reduced from V to V s (e) with 
the definitions 

V m (e) = [k e V : W m {k) > e} and V s (e) = {fc e V : W s (k) > e}, (88) 
where the introduction of the set 7 3 m (e) is actually unnecessary due to the previous remark (II). 



The corresponding phase structure can now be obtained by numerically evaluating equations (|85[) 
and (f5Tj) . For some selected values of the quartic coupling \pj the resulting phase diagrams with 
respect to the parameters Hn and are shown in Fig. [¥] All presented results were obtained for 
an infinite lattice, i.e. L = oo. For —> oo the effective coupling matrix in Eq. (|72[) converges to 
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FIG. 6: Expectation values for the amplitudes of the constant (m: black curve) and staggered (s: gray curve) 
modes for several selected values of the Yukawa coupling constant j/jv and the quartic coupling parameters 
Ajv = 1.0 and Ajv = 10.0. The results were obtained for L = oo. 



the coupling structure of a pure nearest-neighbour sigma-model. One therefore expects a symmetric 
phase centered around kjv = at large values of the Yukawa coupling constant y^ as can be observed 
in the plots. For decreasing y^ the symmetric phase bends towards negative values of k^. In the 
plots the results for the phase transition lines obtained for e = 10 _1 and e = 1CT 3 are compared to 
each other. While the phase transition line to the ferromagnetic phase is unaffected by small changes 
to e as expected, the curves start to differ for the anti-ferromagnetic phase transition at small values 
of j/Af. The discrepancy between these two lines can serve as an indicator down to which value of 
yN the neglection of the modes with W s (k) < e can be considered as a good approximation (besides 
the uncertainties arising from cutting off the power series in Eq. (|5T|) at small values of y^). We 
add here, that we chose the presented parameter range in all phase diagrams such that the volume 
of the space of the considered modes is at least 95% of the volume of the whole mode space, i.e. 
Vol('P s (e)) > 0.95 • Vol("P). For §n — > oo the volume of the neglected modes vanishes and the problem 
encountered during the Gauss-integration in Eq. (|T5[) eventually disappears. 

The order of the phase transitions can again be determined by calculating the expectation values of 
the amplitudes of the constant and staggered modes m and s directly from equations ([55]) and (|57|) . 
The corresponding results are presented in Fig. [5] and Fig. [SI One clearly sees that the occurring phase 
transitions are of second order as one would also expect from the limit j/jv — > oo where the model 
becomes a sigma-model. 



V. SUMMARY AND CONCLUSIONS 



In this paper we have studied analytically the phase structure of a chirally invariant lattice Higgs- 
Yukawa model, originally proposed by Luscher. This was possible in the large A^-limit for small 
as well as for large values of the Yukawa coupling constant and it could be shown that the model 
possesses a rich phase structure. 

In Section [TlTl we began by considering the model at small values of the Yukawa and quartic coupling 
constant and argued that taking only the constant (m) and staggered (s) modes of the Higgs field 
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into account is sufficient for the determination of the phases in that regime of the Yukawa and quartic 
coupling constant. We then presented an explicit expression for the effective potential at tree- level 
in terms of m and s and showed the corresponding phase diagrams for some selected values of the 
quartic coupling constant. In these diagrams all possible phases, i.e. symmetric (m = 0, s = 0), 
ferromagnetic (to ^ 0, s = 0), anti-ferromagnetic (to = 0, s =/= 0), and ferrimagnetic phases (to =/= 0, 
s =/= 0), could be observed. Furthermore, we concluded from our result for the effective potential 
that the occurring phase transitions from the symmetric to the ferromagnetic and anti-ferromagnetic 
phases are of second order. 

In the following Section IIVI we proceeded to the regime of large values of the Yukawa coupling 
constant y^. We showed that for sufficiently large values of un and arbitrary values of the quartic 
coupling constant Xn the model becomes an 0(4)-symmetric, non-linear sigma-model in the large Nf- 
limit up to some finite- volume terms. In particular, this relation to a sigma-model has the consequence 
that a symmetric phase also exists at large values of the Yukawa coupling constant. We determined the 
phase structure of the latter sigma-model by an additional large TV-limit with N denoting the number 
of Higgs field components here. The corresponding phase diagrams revealed again a rich structure 
consisting of symmetric, ferromagnetic, and anti-ferromagnetic phases separated by second order phase 
transitions. The symmetric phase, however, was shown to emerge only in the infinite volume limit. 
For small lattices, finite volume effects cause an asymmetry in to and s which one would not expect 
in a pure sigma-model. These finite volume effects may easily give rise to a misleading interpretation 
that a symmetric phase at strong values of the Yukawa coupling constant does not exist. However, on 
sufficiently large lattices the symmetric phase should become clearly observable and the asymmetry 
should disappear. 

The validity of our analytical results and in particular the latter predictions about the symmetric 
phase at large yjv will be confronted in an upcoming paper with the results of corresponding Monte- 
Carlo simulations including the chiral invariant fermions in a fully dynamical fashion. 



APPENDIX A 



In this appendix we would like to make up for the neglected derivation of Eq. (|56p . Starting from 
Eq. ||53J one finds 



Tr[AB-r = E r ^K»([^- 1 ] Bl|B> -...-[^- 1 ] 



(Al) 



E E 



E 



ipi(ni-n 2 ) ip r {n r .-n 1 ) 



rii,...,n r Cieifci,.-. ] Cre-rfer pi,...,p r £V 
1 e 1 k l ' ■ ■ ■ ■ ** r e r v 



L 4 



L 4 



x Tr. 



8x8 



a^(p 2 ) (S^)($; 3 /|$„ 3 | 2 )) UW,k' 2(p2) V . ... . U Crerk r(pr)a e r{pr) 



Crt n k r £' r e' T k' T 



E E Tr8x8 

ni,...,n r pi,...,p r eV 



n 



a ip i (n i -n i+1 ) 



L 4 



1$, 



, +1 rM( K )(B(^($; i+1 )) u( Pi , Pi+1 ) 



where the definition of the spinor basis transformation matrix E/(pi,P2) given in Eq. (|25[) was used. 
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APPENDIX B 



In this appendix we want to deal with the possibly non-positive eigenvalues of the operator — n e ff+X, 
which would not allow the option of performing the Gauss- integration in Eq. (|75[) over all modes, in 
a more precise manner. We therefore restart our calculation beginning in Eq. (|73|) . Now we perform 
the Gauss-integration solely over those modes O^fc^ir which have their corresponding 

eigenvalue of the operator —n e ff + A not smaller than 2e > 0. We denote the subset of these modes 
as V(e, A). According to Eq. ijTTjl it is given as 



V 



M) = { 



k e V 



-K]y(p 2 COs(fc p 



^2— 2 • q{k) 



A 



Performing the Gauss-integration only over these modes the action reduces to 

.eft ■ nM-^ 2 



> e 



}■ 



(Bl) 



S[m\s l 



. A. a[\ = -ln[dct"(- K e ^ + A)]" JV/z + ^|f][m l ] 2 .(0|- K e ^ + A|0)-L 4 A 

N N ~\ 

EM a -<»l-«^+^k>+E E [4] 2 -(H-- e// + A|fc) 

i=l i=l fce-p(«,x) J 

Tr" In + Al + 2- ■ m 2 ■ L 4 ■ ( -8k N ^ 2 - ^q(0) + x] 

tN V Vn<P J 
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//V 



•iKNip 



vlitp 2 



+ E J-- a l- Li - ( -2K7W 2 ^cos(fc M ) 



16 P 2 / 7 \ X 



(B3) 



where crjL denote the amplitudes of the excluded modes with fc 6 P(e, A), ^ k ^ it and V(e,X) 
V/V(e, A) is the complement of V(e, A). Here the notation 



-N ak 



(B4) 



was introduced correspondingly to Eq. (|79p and the plane wave modes |fc) were explicitly given in 
Eq. (|76p . The determinant det and the trace Tr , respectively, are now only performed over the 
modes k 6 V(e,X), 7^ k 7^ n, as desired. The resulting gap equations can now be obtained by 
differentiating the effective action with respect to to, s, A and all a%. This leads again to Eq. (f8Tj) and 
Eq. (l82|) . Only the third one, Eq. |83|) . is modified yielding now 



2 , 2 

TO + S 



ie?(i,A) 



tep(<,i) 



-KAr^ 2 COs(fc /J ) 



■^2— j«(fc) 



(B5) 



Furthermore, one obtains one additional gap equation for every mode € P(e,A), 7^ fc 7^ tt according 
to 



= a fc 



4 16p 2 

-2k N (p 2 ^ cos ( fc /u) + -2 /n2 



>u=l 



2/^ 



Vfce^(e,A), O^fc^Tr. 



(B6) 



Again we consider the scenario of a purely ferromagnetic phase and the scenario of a purely anti- 
ferromagnetic phase for the investigation of the phase structure. The only particularity here is that 
we assume all to be zero in both cases. (In principle, with this approach one could also study the 
phase structure of some of the amplitudes cr^, but this is beyond our interest here.) We thus arrive 
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directly at the prior equations ([54")) and (J55J, respectively, fixing the value of A as before. With this 
fixation of A the subset V(e, A) now becomes P m (e) for the ferromagnetic phase as already defined in 
Eq. ([55)1 . For the anti- ferromagnetic phase it becomes T' s (e). We have now arrived at the final results 
for the self-consistency equations that were already presented in Eq. [[55]) and Eq. (|57|) . 

In order to get a rough estimate about the validity of neglecting the modes k £ V(e, A) one should 
check the volume of this subset and compare it to the volume of the full set V as we did in our 
discussion in the main text. 
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